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Abstract 

The general properties of a perfect relativistic fluid resulting from the quantum gravitational anomaly 
are investigated. It is found that, in the limit of a weak gravitational field, this fluid possesses a polytropic 
equation of state characterized by two universal constants: the polytropic constant and the natural 
polytropic index. Based on the astrophysical data, the estimates for the polytropic constant are given. It is 
shown that this fluid can describe a considerable part of the cold dark matter. The quantum theory of such 
a fluid is constructed in the framework of the background field method. The Ward identities associated 
with the entropy and vorticity conservation laws are derived. The leading gradient corrections to the 
pressure of the perfect fluid are found and the restrictions on their form are obtained. These restrictions 
guarantee, in particular, the absence of ghosts in the model. The second order nonlinear corrections to 
the equations of motion of a perfect relativistic fluid are analyzed and the explicit expressions for the 
transverse and longitudinal perturbations induced by a sufficiently strong sound wave are obtained. A 
dynamical solution to the problem of time in quantum gravity is proposed. 


1 Introduction 


There are many approaches to explain the general features of the phenomenon of a cold dark matter. The 
common point of view is that the dark matter is a gas of weakly interacting massive particles (see, e.g., 
ii); another one is that the dark matter is a manifestation of deviations of the gravitational laws from 
that dictated by General Relativity (see, e.g., @); the third approach postulates a modification of dynamics 
at small accelerations or introduces the additional tensor fields into the model (see H) ; and so on. In the 
series of papers ||7-I2|, it was proved that the quantum gravitational anomaly exists]^. Its existence leads in a 
natural way to the contribution to the energy-momentum tensor with the properties resembling a cold dark 
matter |^. Such a dark matter is described by the equations of relativistic hydrodynamics with a certain 
equation of state. In this paper, we shall study the general properties of this fluid, with the assumption that 
it is responsible for the most part of a dark matter. We shall also construct its quantum theory using the 


background field method |I4, 151. The latter solves in a dynamical manner the issue of dependence of the 
quantum field theory observables on a choice of the vacuum state for an arbitrary gravitational background. 
In other words, this solves the so-called problem of time in quantum gravity |pb|-p^. 

As known, besides the perturbative nonrenormalizability, quantum gravity possesses another one impor¬ 
tant problem - its precise formulation depends on a choice of the vacuum state for quantum fields. Different 
vacuum states result in different, and in many cases unitary inequivalent, quantum theories. For example, 
it is notorious the dependence of the effective action on a vacuum state in black hole physics 111, 11,1113. 
Usually, the three types of vacuum states are distinguished for quantum fields evolving on the black hole 
backgrounds: the Boulware vacuum | 


I, the Hartle-Hawking vacuum [^, Q, and the Unruh vacuum |R|. 
However, it is clear that one can construct infinitely many vacua leading to infinitely many inequivalent quan¬ 
tum field theories. The same effect, but formulated in other terms, arises as the time problem in quantum 
gravity in-ill- 

A unitary inequivalence of quantum field theories having the same classical limit reveals itself in the 
infrared and ultraviolet divergencies in the effective action. Of course, some divergencies appearing in 
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'^This anomaly is not of the type found in |l^]. 
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the perturbative calculations may cancel upon summing an infinite number of diagrams, but the complete 
cancelation of divergencies in the model without any large symmetry, for instance, in the standard model, is 
improbable. Therefore, one may regard the perturbative divergencies as a measure of a unitary inequivalence. 
The choice of the regularization scheme and concrete values of the constants at the singular structures (by 
introducing the counterterms to the initial classical action) defines a unique representation of the algebra 
of observables. For example, the standard ultraviolet divergencies in quantum field theory correspond to a 
unitary inequivalence of the Heisenberg representation and the interaction picture [^]. For the renormalizable 
theories, a passage from one representation to another is reduced to a renormalization of the coupling 
constants^ (see, e.g., p^ ). 

On including gravity into quantum field theory, the divergencies become dependent on the vector field 

(see, e.g., [^, |^, |^, dehning the physical regularization p4-36| by the energy cutoff, determining the 
Hamiltonian of quantum fields and the vacuum state. If one appeals to the Gelfand-Naimark-Segal (GNS) 
construction then this vector field specifies the vacuum state and, consequently, the representation of the 
algebra of observables (see, e.g.. 


I 37|). Of course, one can choose a certain unique representation fixing 
the vector field by hand as it is done in describing the quantum effects near black holes. However, the 
presence of the non-dynamical field in the effective action leads inevitably to a quantum gravitational 
anomaly, i.e., to a violation of the general covariance (the background independence). One of the possible 
solutions to this problem is to define dynamically. Then the general covariance is restored, but at the cost 
of the additional degrees of freedom. It turns out Q that the dynamics of the field are unambiguously 
determined from the requirement of the fulfillment of the Ward identities for the effective action and are 
described by the hydrodynamic equations of motion. We shall take into account these observations in the 
condition 3 of Sect. concerning the structure of the classical action for the field 

From the renormalization theory viewpoint (see, e.g., |^), one could suppose that the divergencies 

mentioned above ought to be simply canceled out by the appropriate counterterms, and, thereby, the general 


covariance could be restored in the effective action. The principal result of the paper |12| consists in the 
proof that the effective action of quantum gravity contains the finite contributions that are not expressed in 
a covariant form in terms of the metric alone and cannot be canceled by the counterterms. Even the explicit 


expressions for these contributions were found in [1^| for a stationary slowly varying in space gravitational 
background. These contributions are not analytic in the gravitational constant and momenta. The fact that 
such terms should be present in the effective action was frequently mentioned in the literature (see, e.g., 
[O, E^). The non-covariant terms analytic in the gravitational constant and momenta also appears in 


calculating the quantum effects on the black hole backgrounds |27, 41-45|. However, in the latter case, it 
was not clear whether these contributions must be canceled by the counterterms or not. The presence of 
the non-covariant terms in the effective action means that the general covariance is violated in quantum 
gravity. This anomaly is characterized by the vector field referred above, that is, the effective action can 
be made formally covariant, if one introduces this vector field Therefore, following the general 

prescriptions of renormalization theory [ ^ ] and of the effective held theory approach |^, ^|, we have to 
include into the initial classical action of a theory and make this held dynamical. This, in particular, 
implies that a complete cancelation to zero of all the divergencies depending on is unnatural. The various 
examples of quantum held theory models teach us that the structures having a small number of derivatives 
and appearing as divergencies constitute the classical action. 

The fact that the presence of the vector held in the effective action is revealed as the presence of 
some perfect huid resembling a cold dark matter in the weak gravitational held limit was mentioned in 
[p|-p^. In this paper, we develop this proposition. In Sect. |2|, we briehy formulate the main properties of 
the quantum gravitational anomaly. In Sect. the Fock-Taub formalism for relativistic hydrodynamics 
resulting from the gravitational anomaly is expounded. Then, in Sect. Q, we discuss the quantization of 
relativistic hydrodynamics employing the background held method. Actually, we study the dynamics of small 

we establish the general properties of the 


huctuations over the classical background. In Sects. O and 
huid pressure coming from the gravitational anomaly. In particular, we show that this huid possesses the 
polytropic equation of state with a universal polytropic constant and a universal natural polytropic index. 


^Usually, one considers a class of the unitary inequivalent theories possessing certain symmetries; the Poincare invariance, 
the gauge invariance, etc. 
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In the next Sects. [4.3|, 4.4, and [4.5|, the influence of the leading gradient corrections to the fluid pressure on 


the dynamics of small perturbations (the phonons) is investigated and the restrictions on the form of these 
corrections are found. Also, we derive in these sections the explicit form of the Ward identities associated 
with the entropy and vorticity potential conservation laws. Sect. |5| is devoted to the leading nonlinear 
corrections to the dynamics of the longitudinal and transverse phonons. In particular, we find the explicit 
expression for the transverse perturbations induced by a longitudinal sound wave. Then, in Sect. ^ we 
consider a concrete form of the fluid pressure that follows from the renormalization theory. In Sect. 0, we 
discuss some quantum effects induced by the phonons and show that the standard in-out perturbation theory 
has a very narrow range of applicability in the weak field limit 
We shall use the conventions adopted in |1T| 




1/^7 


p _ ryOL 


R = M, 


( 1 ) 


for the curvatures. The square and round brackets at a pair of indices denote antisymmetrization and 
symmetrization without 1/2, respectively. The Greek indices are raised and lowered by the metric which 
has the signature —2. The system of units is chosen such that c = h = 1. 


2 Quantum gravitational anomaly 


In this small section, we summarize the results of the papers |^]^| necessary for us. It was shown in this 


series of papers that, in the background field gauge |14, 15|, the renormalized in-in effective action of the 
helds and has the form (see, for details, [0|) 




( 2 ) 


where <I> denotes all the matter fields including the Faddeev-Popov ghosts both for the matter and metric 
helds. The dependence on the held characterizes the dependence of the effective action on a choice 
of the vacuum state for quantum helds (for instance, using the GNS construction) and, in particular, the 
regularization procedure. 

The requirement of invariance of the effective action with respect to general coordinate transformations 
implies the Ward identity 


The plus sign at the connection means that it is constructed using the metric . The approximate equality 
denotes the fulhlment of the Ward identity on the solutions to the quantum equations of motion of the helds 
$. Henceforth, for brevity, we omit the plus sign at the “-I-” helds. As we see, the Ward identity for the 
matter helds, i.e., a vanishment of any matrix element of the energy-momentum tensor divergence operator, 
is violated. For it is satished, one needs to demand 


= 0 . (4) 

These are the equations of a hydrodynamical type for the held They can be rewritten as 

= 0, eV[^(r,]/u;) = £^ir,/w) = 0, (5) 

where w := The same equations hold for the held too. In other words, the regularized operator 

equality 

V%,[gp,,e,^=0, ( 6 ) 

where all the quantum helds are expressed in terms of the creation-annihilation operators, should be the 
Heisenberg equation of motion of the quantum held 

Despite the fact that, in using the standard procedure to calculate the averages of quantum operators, 
the dependence on the held appears only in the quantum corrections to the effective action, the general 
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scheme of renormalization theory (see, e.g., |34|) says that the field must be included into the initial 
classical action as well. The equations of motion for this field should be equivalent to the condition of the 
covariant divergenceless of the energy-momentum tensor. In this case, a consistent quantization of the fields 
$, and results in the quantum dynamics preserving the Ward identity (Q), (|^. In constructing the 
classical Hamiltonian of the system and in introducing the physical regularization (the energy cutoff), the 
field ought to be used. This field is an average of the operator As expected, we see that the violation 
of the Ward identity leads to new degrees of freedom in the theory - the dynamical field 

The natural normalization conditions are imposed on the effective action involving the field 

1) The Lorentz-invariance. For the flat space-time, where the dependence of the effective action 

on the field should disappear; 

2) The compliance with the Einstein equations. In the flat space-time limit, the vacuum expectation value 
of the energy-momentum tensor operator of the matter fields and the field must be zero; 

3) Minimality. The initial classical action of the field contains only those structures that arise as diver¬ 
gencies in calculating the quantum corrections using the physical regularization (see, e.g., by the 

energy cutoff. It is these terms that specify the energy-time anomaly ||^, and do not include the 
monomials invariant with respect to the global dilatations of 

The last requirement implies, in particular, that, in the flat space-times limit, the initial classical action of 
the field should contain only renormalizable vertices provided the theory of the fields <I> is renormalizable 
in a flat space-time. 

These requirements restrict the form of the classical action of the field but, unfortunately, do not fix 
it unambiguously. When the energy cutoff regularization is employed, the regularization parameter A always 
enters the effective action in the form of the combination ^^/A^. Therefore, if the theory does not contain 
the coupling constants with a negative mass dimension as, for example, the standard model without gravity, 
then it follows from the third requirement above that the coupling constants at the vertices in the action 
of the field possess nonnegative mass dimensions provided the field has been stretched as —)• A^^. 
The inclusion of the gravitational corrections makes this statement improper (see, e.g., p^]). Nevertheless, if 
A <C Mp;, p <C A, where Mpi is the Planck mass and p is a characteristic momentum flowing in the Feynman 
diagrams, then these gravitational corrections are small and can be taken into account systematically using 
the effective field theory techniques |^, 471. In the papers iim, the one-loop divergencies were found. 
Their structure agrees with the property mentioned above. Note that the dilatation —>■ A(^^ allows one to 
get rid of the regularization parameter A in the theory. In that case, it is replaced by the field with the 
mass dimension —1. 

In order to impose the additional restrictions on the form of the effective action and, consequently, in the 
leading order of the perturbation theory, on the form of the classical action of the field one can demand 
the fulfilment of the conformal (Weyl) invariance of the effective action on the cutoff scale A. In virtue of 
the fact that the conformal and the energy-time anomalies are related 10 . this condition fixes uniquely the 
coefficients at the logarithms (see ITlIl and (^^ below) in the effective action of the field : these coefficients 
should coincide with or be proportional to the renormalization group scaling functions /3 and 7. However, 
these renormalization group functions are calculated on the scale A, and all the particles with masses m <C A 
make contributions to the values of these functions. As long as we do not know whether the standard model 
is a final one or there exist more massive particles that are not included into it, these scaling functions are 
unknown. So, this restriction on the coefficients of the action for the field is not constructive in the case 
when the complete Hamiltonian of the theory is hidden. Even if the theory were completely given, it could 
become non-perturbative on the scale A, and the standard methods to determine the scaling function would 
be inapplicable. Therefore, in this paper, imposing additional constraints on the coefficients at the terms in 
the action of the field we shall proceed from the empirical reasonings adjusting these coefficients in such 
a way that they provide the required asymptotics of the energy-momentum tensor in the weak field limit. 

We shall discuss in details the form of the classical potential of the field in Sects. ^ and p. At the 
present stage, it is sufficient for our purposes that this potential obeys the conditions 1 and 2. The explicit 


form of the gradient terms and their influence on the dynamics of the field will be considered in Sect. 4.3 


4 







3 Fock-Taub formalism for relativistic hydrodynamics 


As we have seen in the previous section, we need to construct a model of quantum relativistic hydrodynamics 
in order to obtain a self-consistent quantum field theory involving the quantum gravitational field. At the 
present moment there are many examples of quantization of relativistic hydrodynamics in the literature (see, 
e.g., 1^, ^-^1). However, these approaches are based on the classical action different from that dictated 


by the quantum gravitational anomaly. It is well-known that, upon quantization, different action functionals 
leading to the same equations of motion may result in different quantum theories even for the models with 
finite degrees of freedom (see, e.g., |5£ 


The Fock-Taub formalism |M, is the most appropriate (and 
seems to be unique) to treat relativistic hydrodynamics stemming from the quantum gravitational anomaly: 
the fluid flow is characterized by a single vector field and the variation of the action leads to the equations 
of motion equivalent to the condition of the covariant divergenceless of the energy-momentum tensor (|5|). 
Other Lagrangian descriptions of relativistic hydrodynamics can be found, for example, in 541. 

Let x^(>r), /i = 0,3 and x = {r, a*}, f = 1,3, be a orientation preserving map of the four-dimensional 
manifold V to a four-dimensional space-time M equipped with the metric g^u. Let a = 0,3, be a 

vector field on N such that p°‘habP^ > 0, where hat '■= daX^d^x ^is the induced metric on N. Then the 
action functional |^, (cf. 


S[x{>c)]= d>c^/\h\p[p‘^,hab\, 

Jn 


( 7 ) 


describes perfect isentropic relativistic hydrodynamics in the Fock-Taub representation. Here p is some scalar 
constructed in terms of p^ and hab and their derivatives. As a matter of fact, it also depends on the fields 
<L(x(x)), but this will be irrelevant for our further considerations provided the fields <L(x(x)) obey their 
equations of motion. Note that p°’ is not a dynamical field. It is not varied in taking a variation with respect 
to x(x). 

Since the action depends on x(x) only through the metric h^b and the fields d>, we have 


-—^—= / dx det e(((x ) - )~?—^^ = dete(( 

feM(x) J ^'6gxp{xi>€')) > 5x^^{ 


X 


-29a 


\ 5S ^ 

9/j.p ) + ~ dp9xp ) (8) 


^9\ 


^9Xp 


on the solutions to the equations of motion of the fields d>. The notation eH ■= daX^ has been introduced, 
and we have used the relation 

dbiCpdeteP) = 0. (9) 

The vectors ea constitute a holonomic basis, but not a tetrad. The expression in the square brackets in (^ 
is expressed via the divergence of the energy-momentum tensor defined in the standard manner (^). Hence, 
we arrive at 

= detePy/lilVxT;: = vI^VaT^^ = 0, (10) 

i.e., the equations of motion following from the action ( 0 ) are equivalent to the condition of covariant 
divergenceless of the energy-momentum tensor constructed by the use of this action functional. 

Putting = p^Ca, we obtain the action functional for the field discussed in the previous section. As 
long as 


= e’^Jdbe’^, - Fg^e^) + = 0, 


where F^^ is the Levi-Civita connection constructed in terms of the metric hab, we deduce from (K), 

„ SS 


' 5xP{>c) 


=: vWK = vWVbPa = 0. 


( 11 ) 


( 12 ) 


Performing an infinitesimal general coordinate transform of the variables x in the action (|^, we prove, as 
in the previous section, that 


V^Tab = CpSa + y'^PbSa = d ^ Va(/u;) = 0, Cp{Sa/w)=d, 


(13) 
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where Sa is defined similarly to and the enthalpy density w = p°‘Sa- If one defines the reciprocal 
temperature and entropy density as |P| 


/3 := 0 - := Pw, 


(14) 


then the first equation in the system (1^) is the entropy conservation law for a perfect fluid. It is obtained 
by convolving with p°' the condition of covariant divergenceless of the energy-momentum tensor. 


The second equation in the system (13) implies the vorticity conservation law. The one-form, 

Ta ■■= Sa/w, p°-ra = 1, 


(15) 


is conserved along the fluid flow according to the second equation in (^). For a particular case, when the 
Lagrangian density of the action (|^ does not contain higher derivatives of the fields x^{k) (see, for instance, 
(1^), the one-form (15) is equal to the Tolman temperature one-form ta = Paj |0- The vorticity tensor 
can be defined as (see, e.g., p^) 


Aab := K'^bd[crd] = {dr CpV A r) 


a 6 ? 


Ti'a := - p’^ra- 


(16) 


It is zero if and only if the one-form is integrable. On the solutions to the equations of motion ([T^, the 
vorticity reduces to the external derivative of the one-from and so the latter can be called the vorticity 


potential. Owing to the property of the Lie derivative, the second equation in (1^) entails the vorticity tensor 
conservation. 

In the adapted system of coordinates, where the vector held p^ is straighten [/?“ = (1,0,0,0)], the action 
functional is invariant under the symmetry transformations (see, e.g., 


T ^ t' = T + /(u), a ^ a' = CF'{a)] 

5eP°‘ = 0, SeX^ir, a) = s'"{a)daX^{ t, a), 


(17) 


for any hxed metric gpi,{x). Notice that these local transformations are independent of time r and so are 
not gauge transformations (see the dehnition in ||70| ). Applying this symmetry transform to the action (j^, 
we come to the relation 

'\/\h\Ea = drfa[x{^)\, (18) 

where fa[x{}<)] is a certain local expression. Consequently, we have the conservation laws provided the 


equations of motion are fulhlled. Comparing it with (13) in the adapted system of coordinates, we see that 
the zeroth component of the conservation law is responsible for the entropy conservation, while the spatial 
components are for the conservation. Notice that the zeroth component of the second conservation law in 


(13) is trivial in this system of coordinates. 


3.1 Hamiltonian 

In order to construct quantum relativistic hydrodynamics, we shall not use the canonical formalism explicitly. 
In our case, it is rather awkward since the momenta cannot be expressed through the velocities in an explicit 

J). Quantizing relativistic hydrodynamics. 


form, although it can be done in principle (see for details |31 


we shall use the background held method 0, [TsI which states that the knowledge of quantum evolution of 
small huctuations over an arbitrary background is sufficient to reconstruct the whole quantum dynamics. 
Nevertheless, it is important to realize what the canonical Hamiltonian for the action is. Naively, one 
may expect that it is equal to the energy of a huid. However, we shall see that it is not so for the helds x{>c). 

For simplicity, we restrict our considerations to the case when the Lagrangian density of the action (j^ 
is independent of the higher derivatives (of the second order and higher) of the fields x^{x) with respect to 
r. Then the density of the canonical Hamiltonian takes the form 


d 


Hcan[x{>c)] = x'^ — {y/\hip) - ^/^h[p, 


x^ := a 


(19) 
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Introducing the spatial metric |65-67, 


c[ij .— hij 


x^g^lue'ix^gxpe^j 

x^^g^lux’' 


we obtain 


det hab = x'^ det qij , 


Hcan[x{>c)] = 


( 20 ) 


( 21 ) 


The expression x^dp/dx^ determines a variation of p under the action of an infinitesimal dilatation of the 
variable x^. The local dilatation transform x^ — X{x)x^ acts on the metric hab as 


A(j<) 0 

0 l3x3 


hab ^ ^ab — ^a^cdJb ) — 

As long as p is a scalar, we have 

Pip^^Kb] =P[JbP’'^hab] =p[X{x)p'',hab], 

where the last equality is valid in the adapted system of coordinates. Therefore, 






S[X{>i)p] - 5[p] = J d>cy^\{p[X{x)p°',hab] -p[p°‘,hab]) ~ J dx^/\h\{X{x) - l)i:^^, 

where it is assumed in the approximate equality that A(x) —>• 1. To put it in another way. 


6S 


w 


= = - 


1 . Bp 


^x^^ ’ 

in the adapted system of coordinates. As a result, the density of the canonical Hamiltonian reads 

Hcan[x{x)] = -x/\q\ldw = -a/Mct, 


( 22 ) 


(23) 


(24) 


(25) 


(26) 


that is, in the Fock-Taub formalism, the canonical Hamiltonian for the helds x{x) coincides with the total 
entropy taken with an opposite sign. This result can be foreseen from the following considerations. The 
equations of motion (|T^) possesses a conserved charge coinciding with the total entropy of a fluid. Inasmuch 
as the Lagrangian ( 0 ) does not depend on r explicitly, the corresponding canonical Hamiltonian is an integral 
of motion too. Hence, in a general position, this integral should be proportional to the total entropy. One 
can also add to the integral of motion a linear combination with constant coefficients of the integrals of 
motion following from (^). However, the only natural constant vector in the model at hand is (in the 
adapted system of coordinates) and we saw that p°‘fa is proportional to the entropy density. 


Formula (^) is a generalization of the well-known thermodynamic relation for the canonical ensemble 


odp 


(27) 


where p is the pressure. Carrying out the standard reasonings analogous to those used for the systems with 
the Hamiltonians equal to the total energy, one can see that the part of a fluid characterized by bounded 
values of the coordinates cr* on the manifold N tends to increase its entropy, i.e., to diminish the magnitude 
of the canonical Hamiltonian corresponding to this part of fluid. The phonons with the mode functions 
diagonalizing the quadratic part of the Hamiltonian take away the negative entropy to spatial inhnity. This 
leads to increasing of the entropy of the given part of a fluid. The mention should be made once again that 
the total entropy of a whole fluid is conserved. Other more effective mechanisms to produce the entropy in 
a perfect fluid are also possible. For example, it is produced with the help of shock waves or a turbulence 
(see for details, e.g., 0), when the map x{x) ceases to be one-to-one. However, this issue needs a further 
investigation. 

The density of a canonical Hamiltonian in the space-time M, i.e., for the helds x“(x), coincides with the 
energy density as it should be. This canonical Hamiltonian is not conserved for the case of the non-stationary 
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matter and gravitational fields. Indeed, in virtue of the general covariance, the action ( 0 ) can be rewritten 
in the form 

SMx)] = ^dxA/|9|(x)p[/(x(x))e^,5-^^(x)], , (28) 


i.e., as a functional of the fields x:(x) defined on the manifold M. The equations of motion reads 


= = - vl^v.r; = o. 


(29) 


where is defined similarly to T^. Notice that, in taking a variation with respect to >c{x), the argument 


of is also varied. If the action functional (28) does not contain higher derivatives of the fields x(x) then 
employing the relation, 

det 


= -ea4, 


(30) 


it is easy to show that the density of the canonical Hamiltonian, 

'-{x)] = -p) = +p)= \/\g\{wu^UQ-p) = (31) 


Hr. 


coincides with the standard expression for the energy density. Here is the four-velocity of a fluid 

flow. Henceforward, we shall consider the Fock-Taub representation of relativistic hydrodynamics in terms 
of the fields x(x) rather than x{x) since the symmetry transform dT?] ) looks simpler in this representation. 
However, it may happen that the Fock-Taub formalism in terms of >c{x) is more appropriate for some 
problems. 


4 Quantization of hydrodynamics 


As we have already noted, according to the background field method |14, 15|, to construct the quantum 
field theory, it is sufficient to know the quantum dynamics of small fluctuations on an arbitrary background. 
The quantum theory of a non-relativistic fluid is also developed in this way |f7^ ]. Thereupon our task will 
be to describe the dynamics of small fluctuations (the phonons) of a relativistic fluid taking into account 
the conditions 1-3 given in Sect. The one-loop calculations performed in [^, |ll], Q show (see the 

condition 3) that the action of hydrodynamics (^ is highly nonlinear. Therefore, we expect that quantum 
hydrodynamics is perturbatively nonrenormalizable. In order to conduct consistent calculations, we stick to 
the effective field theory approach (see, e.g., |4f:, ^) supposing that the gradients of the background fields 
are small, i.e., in fact, considering the infrared limit of a theory. 


4.1 Classical potential 

If we neglect the gradients of all the fields and take into account that the only known macroscopic fields with 
nonzero vacuum expectation values are the metric the Higgs field </>, and the vector field (we also 
neglect the small, but possibly non-vanishing, contributions depending on the electromagnetic field) then, 
from the condition 3, the expression for the fluid pressure follows (see, e.g., 0 ) 


v«(p\ = E E In* 


(32) 


fc = 01=0 


where ta ■= Pa/p^- Such a peculiar structure of the classical potential is caused by the fact that the 
regularization parameter of the energy cutoff regularization always enters the nonrenormalized effective 
action in the combination and the general structure of divergencies has the form (32) provided the 

terms suppressed by the Planck mass are omitted (see the remark in Sect. § after the condition 3). The 
quantities a^i are the curvature independent gauge invariant scalars with a mass dimension 4, 


2-fc 


aki{\4>\‘^) = ^ bkis\4>\'^'', bkis = const, 


(33) 


s=0 
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that is, is a dimensionless quantity. The logarithmic corrections describe the so-called anomalous scaling 
(see, e.g., |73|) on the cutoff scale A and are assumed to be small. The whole expression 


specifies the 

energy-time anomaly 00 of the effective action so long as the gradients of the fields are discarded. 

The flat space-time limit corresponds to = 1. Indeed, for a weak slowly varying gravitational 

field, the vector field is close to the Killing vector of the corresponding stationary metric approximating 
the slowly varying one since in a stationary case the Killing vector field is a solution to the equations of 
motion ( 0 ) i |4^ . So, in that case, we have in the adapted system of coordinates (see, e.g., 

f = p‘^ = 


1 -|- 2(/?Ar, 


(34) 


where is the Newtonian potential. Of course, one can normalize the vector field to any other positive 
constant in a flat space-time. In a certain sense, this corresponds to a shift of the Newtonian potential by 
a constant. Furthermore, it is natural to redehne the field as —>■ where A^^^ is some constant 

parameter with a unit mass dimension, such that p“ has the mass dimension 1. Then the entropy density 
of a fluid (14) will possess the standard mass dimension 3. It has the dimension of an energy density for 
the normalization condition = 1 in a flat space-time. However, from an utilitarian point of view, these 
dilatations of the vector p“ just results in the redefinition of the constants Uki and does not change the 
observables. Such transformations resemble the renormalization group transformations and coincide with 


them under the condition that the coefficients at the logarithms in (^) coincide with the scaling functions 
j3 and 7 . Thus we put p^ = 1 in a flat space-time. 


The classical pressure (^) should satisfy the conditions 1 and 2 in the leading order of the perturbation 
theory. Since 

= -2p'p“ 



(35) 

= p“/then 


H|(/>|=r),p2=l ~ 

(36) 


Po{p^ = 0, Po(p^| 

where p ~ 247 GeV is the vacuum expectation value of the Higgs held in a hat space-time. Is also supposed 
that the term ooo in Po contains the classical potential of the Higgs held in a hat space-time limit. Therefore, 
the following relations should hold 


Po{p^ 


')||^|=,,,p2=i = 0, Aif po{p\ = -ml 


(37) 


where the dot denotes the derivative with respect to |(/>p and m 

a hat space-time. The normalization conditions (^), (^^ are only necessary ones. The experimental data 
may lead to additional constraints on po- For instance, it was shown in ||9| that, in order to reach an agree 
of the theory with the experimental data concerning the gravitational red-shift effect, one has to demand 


125 GeV is the Higgs boson mass in 


P'^{p'M?)\\p\=r,,p^=l = ^- 


(38) 


Let us stress once again that the normalization conditions (0), (0), and (^) are imposed on the effective 
action and are valid for the classical potential only in the leading order of the perturbation theory. Having 
calculated the quantum corrections to the classical potential, we need to impose this conditions to hx the 
arbitrary constants appearing in the effective action. It is the effective action which is observed in an 
experiment. 


4.2 Polytropic equation of state 

In this paper, we identify the fluid resulting from the existence of the quantum gravitational anomaly with a 
dark matter or its considerable part. In the limit of a weak slowly varying gravitational field, it follows from 
(36) that wq ^ Eq K, aQ ^ X and pQ ^ x'^, where, for short, we have introduced the notation x := 1 — p^ and 
eo is the energy density of a fluid. In a general case, we can adjust the constants in po in such a way that 
the first n derivatives of po with respect to p^ at = I vanish. Then we have 


Eqk, wq = — 2/3^Pq « 2Ax'^' 


„n+l 


Po 


A- 


n -|- I 


n G IKI, 


(39) 
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where A is a positive constant with the mass dimension 4. Rearranging these expressions, we arrive at 


l+l/n 


PO = 


2 (n + 1 )( 2 ^)V>^ 


=: Ke, 


l+l/n 


eox 


2 (n + l)' 


(40) 


""" den 


This is the equation of a polytrope (see, e.g., [Q) with the index n and the polytropic constant K. If we 
identify Kdmx/[i{n + 1)] with kT then the equation of state ( ^0|) coincides with the equation of state for a 
perfect gas. The limit n —?• oo corresponds to the isothermal polytropic process. The speed of sound in the 
weak field limit becomes ||9-H| 

dpo Po _ ^ 

m 3p(, + 2/32p" ^ 2n' ^ ^ 

Thus, in the weak field limit, the relativistic fluid stemming from the existence of the quantum gravitational 
anomaly behaves as a perfect gas in a polytropic process with a universal polytropic constant K and a 
natural polytropic index n. This fact can be used to falsify the hypothesis about the identification of this 
fluid with a dark matter. 

Notice that, in a stationary case, when the fluid is in a hydrostatic equilibrium, the asymptotics (^), (^^ 
imply that this fluid is absent in a flat space-time limit. This guarantees the absence of a Lorentz-invariance 
violation as we have demanded imposing the normalization conditions. On the other hand, as long as the 
fluid itself produces the gravitational field, there may exist stable self-sustained configurations consisting of 
this fluid and the gravitational field only. The cosmological constant and the cosmological expansion caused 
by it decrease the stability of such gravitationally bounded systems (^, |^. So these stable self-sustained 
configurations may exist only if their density is larger than a certain critical one ||75| , 

We can estimate the value of the constant A and, consequently, of the constant K assuming that the 
velocity of motion of the dark matter with respect to the Galaxy center in the Sun neighbourhood is non- 
relativistic. In this case, we can set x ~ —2(/?q and take |7^^| (see, however, 1^) 


|(^o| ~ 10 


-6 


£0 ~ 0.5 GeV /civ? 


(42) 


whence 


A = 


e© 


2 |2(^o|-- 

The concrete values of A for several different n are as follows 


Ai Rs 1.25 X 10^ GeV/cm^ 
Ag « 9.8 X lO^^GeV/cm^ 


(leV)^ 

(9.3 GeV)^ 


(r//26)^ 


A^k. l.Q X 10^^ GeV /cTc? 
^9 Ri 4.9 X 10^° GeV/cm^ 


(19keV)^ 
(248 GeV)^ 


T] 


(43) 


(44) 


At first sight, it appears that the natural value of n for the electroweak scale is n = 9. However, as we 
shall see, such a value is n = 8 since a small factor of the order of 1 /n! diminishing the magnitude of the 
constant A arises in developing the pressure po a Taylor series in x. For n < 8 , the constants in the 
coefficients a^^d^p) have to be adjusted in such a way as to diminish the energy density and, consequently, 
the magnitude of the constant A from the natural electroweak scale to an acceptable value. 

The models of a dark matter with a polytropic equation of state have been already considered in the 
literature. As well-known, the conventional Navarro-Frenk-White (NFW) prohle for the dark matter 
density in galaxies possesses a cusp at the galaxy center which is not observed in the astrophysical data. The 
observed prohles are more likely described by the so-called cored distributions (for their scaling properties, 
see, e.g., [p^]). The NFW profile fits the results of the N body simulations and the cusp in it is due to a 
pressureless of the dark matter constituted by weakly interacting massive particles (the dust). The polytropic 
equation of state introduces a non-vanishing pressure and so avoids the cusp problem. In the papers |p5|-|^, 
a dark matter with self-interaction was studied. Its equation of state was taken to be (^) with an arbitrary 

, The polytropic equation 


n and the interval 3.5 < n < 5 was approved as the most probable one 


of state also arose in employing the Tsallis formalism to describe the dark matter. The latter was supposed 
to consist of collisionless particles in a non-equilibrium thermodynamical state (see | 


TO I and the critique 


in |91|). These papers suggested the optimal value of the polytropic index to be n ~ 5. In the paper |75|, it 
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was argued that, out of the cusp, the NFW profile for the dark matter density can be obtained by using the 
polytropic fluid with the index n = 3. In describing a dark matter as a Bose-Einstein condensate of some 
scalar field (see, e.g., [p^95|), the polytropic index is found to be n = 1 in the non-relativistic limit. 

The cosmological data can also impose certain constraints on the parameters of the equation of state 
(40). We shall not give here a detailed analysis of the cosmological model with a dark matter of the form 
(ii, &, but only check that the fluid with a polytropic equation of state does provide a model of a cold 
dark matter. Let us pass to the system of coordinates where the fluid at issue is at rest in average and its 
density is homogeneous, i.e.. 


— (^°(t)^ 0,0,0) and = dt^ — a^{t)^ijdx^dx^, (45) 

where ds^ is the interval squared of the Friedmann-Lemaitre metric. At the present stage of the Universe 
evolution, the average density of the dark matter is small [||, ^ 


~ 10 °GeV/cm'^, Sc = 0.53 X 10 °GeV/cm' 


(46) 


where Ec is the critical energy density and ^dm = 0.20 is the fraction of a dark matter in the energy content 
of the Universe. Therefore, we can consider that the fluid is in a weak field regime ~ 1 and obeys the 
equation of state (^). Then the condition of covariant divergenceless of the energy-momentum tensor (1^) 
possesses the solution 


^ — ^dm 


ao 


) {i+AdL” 


1 - 


3/ni —n 


(?)'”]} 


2Ax^^ 


(47) 


where oq is the present-day value of the scale factor. This solution corresponds to a cold dark matter (a 
dust) so long as Ke\^ <C 1. From (|40|) and (|4^), we obtain 


K = 


(2A) 


— l/n 


l<7>©|e 


— l/n 
© 


2 (ra + l) 

Taking into account the magnitudes (^ 
that <C 10“® for any natural number n. 


n + 1 


^ l/n ^ 

n -|- 1 V £0 / 


(48) 


(P^) of the quantities entering the last expression, it is evident 


The solution (L7) allows us to find the cosmological value of q{t), i.e., the value of averaged over the 


homogeneity cell of the Universe. Neglecting in (| 

I i/n/«o\^/"' 

= 2(n + l)K£j^ 


we come to 


\2a) 


CLq N 3/*^ 

a ) 


CL / V £0 / 


l/n 


a ) 


'i/n 


(49) 


One can naively surmise that In^^ = ln(a/ao) (compare, e.g., 1^, P, with @-Q), but now we see 

this is not the case. The solutions of the Friedmann equations for a dark matter with a polytropic equation 
of state for an arbitrary n > 0 can be found, for instance, in 


llOOl, [I01| . The polytropic gas with the 
negative index n was used, for example, in |^, |100| , |102|] to describe a dark energy. It should be stressed in 
this connection that, for the model ( ^^ we consider, the polytropic equation of state (^^, ( ^0|) arises only 
in the weak field limit. Formulas (39), (40) are not applicable at the early stages of the Universe evolution. 

Notice also another one feature of the equation of state (^), (|^). For n odd, the energy density of a fluid 
can be negative, if there are domains in space where x < 0, i.e., > 1. For a spherically-symmetric case, 

the problem of a hydrostatic equilibrium of a polytropic fluid is described by the Lane-Emden equation (see, 
e.g., I?^ ), and it is well-known that such regions exist with a necessity at n < 5. Usually, the transition from 
the region with x > 0 to the region with x < 0, in increasing the radius r, is interpreted as the presence of a 
sharp boundary of a gas cloud. The part of the solution to the Lane-Emden equation with r > xq, where xq 
is the transition point, is thrown away as unphysical: the particle number density is negative at x > xq. This 
procedure is quite reasonable for a usual gas of particles, but, in our case, it would be unnatural to impose 
the additional constraint £ > 0. Usually, the domains with a negative energy density appear in the solutions 
to the Lane-Emden equation at sufficiently large values of x. One should also take into account that, at 
large x, the negative contribution of the cosmological constant to a gravitating “mass” becomes considerable 
since M\^{r) = —87rAcX^/3, where Ac is the cosmological constant. In particular, for the Local 


77 


group of galaxies, the contribution of the cosmological constant dominates at x > 1.4 Mpc |76, 103|. One has 


to take into account the cosmological expansion and the galaxy formation processes on these scales |104|. 
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4.3 Gradient corrections 


In what follows, we shall consider the general case (32) with the asymptotics (|3^ , (^). In the stationary 
limit, the gradient corrections to the pressure (^) were found in |11| to the one-loop approximation. In this 
paper, we restrict our study to the leading gradient corrections that are of the second order in derivatives 
and the homogeneity degree —2 under the dilatations p°‘ \p°‘. One could hope that some structures do 

not appear in calculating the quantum corrections (this does happen for higher coefficients of the heat kernel 


expansion, see [1051). However, all the possible structures permitted by the two restrictions above arise in 
our case. 


Indeed, in jllj, the quadratic divergencies of the second order in derivatives were found in the one-loop 
effective action of quantum gravity induced by a scalar field. Up to the total derivatives, they are 




RV 


where fab ■= d[atb]- The one-loop contribution of fermions produces the additional structure [ 106 , 




(50) 


(51) 


where := da In These expressions were obtained under the assumption of stationarity of the space- 

time, that is, the vector is the Killing vector for the metric hab here. In a general case, we have eight 
independent structures submitting to the restrictions pointed above such that 


P2 = biRt^ + + bsVatbVH^’ + b^iVat^ + + beiaaty + hp^^Vabt^ + (52) 

All the constants hi have the mass dimension 2. One can verify that, in the case of a stationary metric 
with the Killing vector p“, all these structures reduce to the three structures presented above up to the total 


derivatives. Obviously, we can choose other basis of independent structures. However, the basis (^) seems 
to be the most convenient for a subsequent derivation and analysis of the equations of motion. 

Our task is to find the constraints on the coefficients hi that provide 

i) The absence of unstable modes (the ghosts) in the weak field limit; 

ii) The absence of small perturbations (the phonons) in the flat space-time limit p^ = 1, i.e., in this limit, 
they must possess the dispersion law of the form a; = 0 or a; —>■ oo, where oj is the energy of the phonon 
mode. 

The necessity of the first condition is evident. The second requirement is needed for the Lorentz-invariance 
to be preserved in the flat space-time limit, i.e., the condition 1 of Sect. Indeed, so long as the phonons 
are bosons and can constitute a condensate, the existence of nontrivial phonon branches of the dispersion 
law inevitably leads to the possibility of macroscopic perturbations of the field in the flat space-time limit. 
Consequently, it leads to a violation of the Lorentz-invariance. 

Note that the analysis we are about to carry out differs from that is given in the papers |108- 110 |. In 


those papers, the general equations of motion of a relativistic fluid with viscosity were considered and the 
conditions were found when the divergence of the entropy current is nonnegative on the solutions to the 
equations of motion. In our case, the equations of relativistic hydrodynamics are Lagrangian, describe a 
perfect fluid, and the divergence of the entropy current is zero (|T^). 

In order to obtain the contributions of the terms ( |5^ to the equations of motion (1^), we need to find 
the variation of the action ( 0 ) with the pressure (^^. It is not difficult to verify that 


<5p“ = 0, bhab = V{a^b)^ 
5Kb = hi^{ab)0^-R\ab)kO^) 
= -2kt‘^H'^p^Va9b, 


SVWl = vWVa0“, 

= Vbae^ + R^akbO^ 


baa = 


(53) 


where 0“ := e^5x^. Therefore, the variation with respect to is equivalent to the evaluation of the Lie 
derivative along the vector field 0“ with the only difference that the variation of p°' vanishes. Formally, one 
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may put Cqp'^ = 0. Taking the variation of the terms (^) in the action (|^, we deduce 

El = 2Va{t^tcR) - RVcR, El = Vdi^RabtRh^Pc - 2R,atR^) - RabVcitR^), 

El = 2[VhaVcR + - 2tR’^p,V\)] , El = Vc[{VaR? + - AVS’^Rpc^abt'^), 

El = 2[vHV^t + Va{Ru,Vh)], ^2V4tV[Vfe(tVaa) + 2(aat“)2]} - V,[(aat“)2], 

El = 4S/d[th\Vab{Rt^)\ - Vc{p"p^Vabt^) - Vab{p"p^VA " p^^p^Vdt^R^cb, 

El = Ve[pVfe(p“Vat^)] - AVd[th%Va{p^ybP% 

(54) 

where t := \/^ and the index at E^ marks the number of the term in (^). The term without derivatives 

El = -VcPo - 2Va(p“pcPo)- (55) 


(p^ gives 


4.4 Linearized equations of motion 

The background field method is realized in the standard way. The field x{>c) is represented in the form 

x^(x) =: x^(x) + 0'^(x) =: (56) 

where x(>c) is an arbitrary background field that afterward will be identified with the vacuum expectation 
value of the operator x{>c), and Ca = daX^- It is assumed that the vacuum state of a fluid is the ground 
state of the Hamilton operator associated with (|^). Since the Hamilton operator does not depend on the 
time r explicitly, its ground state evolves in a trivial manner when the map x(x) is one-to-one. So one can 
use the in-out effective action to calculate the vacuum averages of operators. From the quantum field theory 
formalism viewpoint, the entropy can be generated by the non-perturbative processes only when the map 
x{x) develops the fold singularities. 

The symmetry transform (^) acts as |]4, ^ 


5,x^^{>^) = 0, 40 “ = e“ + , (V,)^ := dadl + 


(57) 


where Va is a trivial symmetric connection constructed with the help of the Jacobi matrices Ca- The kernel 
of the symmetry generator looks as follows 


40^ = e^R^, 


R^ = {5l + Vj’^{^'ma-a'), 


where a := 


{a, a) and B := {b,>c'). Obviously, these generators form an algebra 


[4i) 42] *^[ei,£ 2 ]‘ 


(58) 


(59) 


The effective action of the quantum fields 0“ depends parametrically on the background fields x{>c). It is also 


invariant with respect to the transformations (IT). This symmetry of the effective action can be violated 
only by the quantum anomalies. However, the physical regularization, that we imply, is covariant under 
the transformations (^) and cannot destroy this symmetry. Also there is not any additional structures in 
the effective action as, for example, in the case of a scalar field on a curved background (see, for details, 
il, [TT| , [l^ , ^3|-^5||). The fluid Hamilton operator defining the dynamics and the vacuum state is 

expressed in terms of the fields and x(x) which are already present in the effective action. Therefore, we 
have the Ward identities 


Ri 


6r 


= 0 , 


Ri 


4 r 




6R^ 5T 


= 0 , 


Ri 


4 r 




SR^ 4 r 


§0A ’ « 69^59^ 69^ 69^ ’ “ 69^69^69^ 69(^ 59^59^) 

In particular, the relation for the polarization operator follows from the second identity 


= 0 , etc. 


fiTn“^(r, cr; 44) = 0 , 


H 


AB 


:= {D 


-l\AB 


4 r 


69^59^ 


6»'4=o’ 


(60) 


(61) 
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on the solutions to the quantum equations of motion. Here >c') is the propagator of the fields 0“ on the 

background x(x). In fact, the Ward identities found express the entropy and vorticity potential conservation 
laws. 

In order to investigate the dynamics of small perturbations of the field x(x), which are determined by the 
classical equations of motion, it is necessary to linearize the expressions (^) and ( |55[ ) over the background 
x{>c) (henceforward the bar is omitted). To this aim, it is useful to employ formulas ( |53|) once again. The 
expressions for the linearized equations of motion are rather huge and given in Appendix We only write 
here the linearized expression (^): 

JEO = -2[Ve(pVVoVa0fe) + Va(pVVo)V(,0fe) + + 0^Rsbca) + 

+ p’^PcP'o^abO^ + 2Va{p’^PcP^P^Po^bes)] • (62) 

This expression will be needed for us in deriving the quadratic in 0“ terms in the fluid equations of motion. 

In the limit when the gradients of the metric field can be neglected and 


it follows from formulas 


/ = /3<5o“, 

and (|125|) that 




= -2/3Vo dcOo + 0c + + 2/3^ 

Po 


,2P0-A 


(4) 


6Et = 2/3-2 - 2d Jo - 2<5O(aa0 “ - 2 0 o)]: 

5El = 8/3-2(5e0O + 


/3 = 'JfP‘ = const, 


6El = 0 , 6E^ = 0 , 


6E^ = 2/3-2,5°n0o, SE^ = 2/3-^5^ Oo, 
5El = - 4 / 3 - 2 ( 9^00 + 


(63) 


6 El = 2 / 3 - 200 ;, 


- 2 e 0 


(4) 


(64) 


where □ := E^^dadh is the d’Alembert operator and the dots denote the derivatives with respect to r. The 
hrst and second expressions vanish not only in the limit referred above, but for any vacuum solution to the 
Einstein equations Rab = 0. The last two expression are similar in the limit (|6^ . 

Let us consider, at first, the linearized equations of motion (|6^ in the limit ([6^). In this case, 

6El = -2/3Jo{dJ + Bi) = 0, SEl = 2jp'o(c7% + W = 0. (65) 

It is convenient to split the perturbation 0j into the longitudinal and transverse components 

Bi = b} + B^, 0 | = diB, diBl = 0 . ( 66 ) 

Then the longitudinal modes obey the standard wave equation describing the propagation of phonons 

0o = -0, 6El = -2pyj-J-M) = ^, (67) 

where we have assumed the zero boundary conditions at spatial infinity. The dispersion law takes the form 
a;2 = Cgk^ and reduces to ^2 = 0 in the flat space-time limit. The dynamics of transverse modes are described 
by 


00 = 0 , 


ii-L 


= 0 , 


SE^; = 2jp'oc;% = 0 . 


( 68 ) 


This is the so-called entropy-vortex perturbations (see, e.g, |71|). They possess a trivial dispersion law 
a;2 = 0, that is, they do not depend on the Lagrangian time r and move along with the fluid flow. The 
existence of such peculiar perturbations is a consequence of the vorticity conservation law (18). 

These entropy-vortex perturbations represent a certain problem for a straightforward application of the 
standard quantization procedure E 


1111 - the propagator does not decrease for large |k| that results 
in severe ultraviolet divergencies, and it possesses a singularity at w = 0 that leads to the presence of infrared 
divergencies. The physical solution of this problem, in contrast to the formal one []T1|, consists in that one 
needs to take into account the gradient corrections to the pressure and nonlinear terms in the equations of 
motion. These corrections modify the dispersion law. 
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The gradient terms of the second order in derivatives ( p^ do change the dispersion law for the transverse 
modes. However, we shall see now that this modification is prohibited by the normalization condition (ii). 
Multiplying the contributions (34) by the corresponding coefficients bi and add them, we derive 


(4) 

2/3 ^ (63 + 64 + 65 + 65 + 867 — 4bs ) 6 q — (63 + b ^) A9Q + (64 — dfey + 2bs)di Oi + / 3 "^Po(c^ ^^0 + di9i ) 


2(3 


-2 


(4) 


bsi 9i - A9i) - b4dij9j - (64 - 467 + 368)5* 0o - /3Vo(5i6'o + Oi) 


= 0 . 


For the transverse modes of the form 

we arrive at 

Whence the dispersion law follows 


9a = ( 0 , 0 *^), 

b^U9i - P%9i = Q. 


= 0, or = —/3 ^Po/^3 + 


= 0 , 

(69) 

(70) 

(71) 

(72) 


In the latter case we have the dispersion law of a relativistic massive particle with a mass squared —/ 3 ^Pq/ 63 . 
Bear in mind the analogy with the sound waves in an isotropic elastic medium (see, e.g., |112|| ), one may 


expect the existence of such a dispersion law for transverse modes when the gradient corrections (52) are 
taken into account. Similar gradient corrections to the free energy of a solid result in the nontrivial dispersion 
law for the transverse perturbations. In the flat space-time limit, this branch of the dispersion law describes 
the massless perturbations and does not die out. This contradicts the condition (ii). Hence, we set 


^3 = 0 , 


(73) 


and the issue of the transverse modes is left. We shall return to this problem in the next section regarding 
the nonlinear contributions to the equations of motion. 

4.5 Stability criteria 


Let ea(cu,k) be the eigenvalues of the Fourier transform of the matrix operator (69) acting on 9°“. Then the 
dispersion law determined by the equation 


ea(wa(k),k) = 0, Vk, 

describes stable perturbations as long as im 


a;a(k) G K, oje'^{uj) 
It is useful to rewrite the latter condition as 


t*;=a;a{k) 


> 0, Vk. 


d£a 


duj^ 


UJ=UJa{k) 


> 0 . 


(74) 


(75) 


(76) 


It ensues the correct sign of the propagator residue at the pole in the a; plane: it is positive for a; > 0 and 
negative for w < 0. The reality condition for a;a(k) can be weaken to some extent to allow for unstable 
particles. Then it is required that 


Rea;a(k) Imu;a(k) < 0 and Reu;Ree(j(a;) 


ljJ=UJaO<-) 


> 0 . 


(77) 


However, this case is not realized for the wave operators corresponding to the Lagrangian equations of motion, 
i.e., to the equations of motion coming from a variational principle. 
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In the framework of the effective field theories, i.e., the derivative expansion that we are employing now, 
the terms with higher derivatives with respect to time (of the third order and higher) in the equations of 
motion cannot produce extra branches in the dispersion law provided the conditions ( |75D are fulfilled. These 
terms just modify, at large energies and momenta, the dispersion law given by the equations of motion of the 
second order in derivatives with respect to time (|65|). Indeed, expanding the action in derivatives, we obtain 
from it the wave operator polynomial in u) and k. This operator possesses the bounded eigenvalues ea(cu,k) 
at finite uj and k. Equation (74) has two roots, at least, for small oj and k. These two roots satisfy (75) and, 
up to small corrections, are determined by the theory without higher derivatives (|65l). Under sufficiently 
weak restrictions on the wave operator of the Lagrangian equations of motion, the eigenvalues £a{oj) obey 
the Schwarz symmetry principle (see, e.g., uni) 


=ea(w*). 


(78) 


So, if (74) possesses a complex solution then a complex conjugate to it is also a solution to (74). Hence, 
for a stable model, equation (|7^ has only two nondegenerate roots (referred above) in the complex oo plane. 
Otherwise, as long as (74) has some additional complex roots, they necessarily violate the condition ([77[). If 
these extra roots are real then the continuity of £a(^) implies with a necessity that the condition ( [75|) does not 
hold for one of these roots, at least. For gauge theories, these considerations are valid for the physical modes 
only. The nonphysical modes can violate the conditions (^), (^). Only when £a{oj) possesses singularities 
on the real axis in the ui plane may the additional stable branches in the dispersion law appear. However, 
these singularities cannot be obtained at any hnite order in the derivative expansion. 


In our case, we need to analyze the system of equations (TO) for the longitudinal modes 6- at 63 = 0. 
Assuming the zero boundary conditions at spatial inhnity, the second equation in (| 6 ^) gives 


e = -(/ 3 Vo + hA)-^ [/SVo^o + (64 - 467 + 2bs)eo] ■ 

Using this expression, we have for the first equation in 


(79) 


2/3 


-2 J 4 
^ id 


64 + 1)5 + 60 + 867 — 46 g — 


(64 — 467 + 21 ) 8 )^ 2 

/3V'o - 


— UJ 


c ^/3Vo + + 2 


/3Vok (64 - 467 + 2bs) 


/3Vo - biid 


+ 


P^p'o - bdd 


9o = 0, (80) 


where 8q = 0o(^)k) is the Fourier transform of Equation (|^ ) has four roots in the complex uj plane 

However, as we have already established, this is impossible for stable models. In order that equation 
has only two nondegenerate roots, we have to demand 


64 + 1)5 + 6 g + 867 — 4^8 — 0 , 64 — 467 + 2bg — 0 . 


With the account of (H), the condition (^) is satisfied as long as 

65 < 0. 

The dispersion law takes the form 


a ;2 = 


(64k2 - /3Vo)(-65c2k2 - ’ 


whence it follows that 


64 > 0. 


(81) 


(82) 


(83) 


(84) 


In the flat space-time limit, the dispersion law (83) reduces to = 0 as it should be. At the small momenta 
k, the terms with higher derivatives are negligible, and we have the standard linear dispersion law for phonons 
with the speed of sound Cg- At large momenta, (x)(k) tends to zero provided that 64 7 ^ 0 and 65 7 ^ 0. Such a 
behaviour of the dispersion law is unphysical since, in that case, there are perturbations with a tiny energy 
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(and a large momentum) possessing a negative inertial mass. We identify our fluid with the dark matter 
that does not display such a property. Therefore, either we set 64 and/or 65 to zero, or suppose that, at 
large momenta, where w'(|k|) becomes negative, it is necessary to take into account the corrections to the 
pressure of the higher order in derivatives than in (^). These terms may correct the unphysical behaviour 
of the dispersion law at large momenta. 

The general solution to the equations (^) and (^l|) can be written as 
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where a are arbitrary constants. Evidently, 


the inequalities (^) and (M ) are equivalent to 


(85) 


C 3 < 0 , C 4 < 0 . 


( 86 ) 


As a result, 

P 2 = 2ciGabt°'t’’ + C 2 Rabt°‘t^ + C 3 pr“^ - 2Ci{V- C 5 Va(/of“Vbp^')t^, (87) 

where 

Gab •— Rab 2^abRj P^ab '— hab UaUb- 

The case 64 = 0 corresponds to C 4 = 0, and the case 65 = 0 implies C 3 = 0. The inequalities (^) assure 
that the structures standing at these coefficients give nonnegative contributions to the pressure. Further 
constraints on the coefficients q can be found, if one analyzes the stability of small perturbations on a 
curved background and nonuniform fluid flow. In this paper, we shall not conduct such a study. We should 
only note that the structures at Ci and C 2 are of a fixed sign and, consequently, can be made nonnegative as 
long as the weak (for the structure at ci) and strong (for the structure at C 2 ) energy conditions are fulfilled. 
The structure at C 5 has not a definite sign, and so one may surmise that a further investigation of stability 
will lead to C 5 = 0 . 


5 Nonlinear corrections 


In this section, we consider the quadratic in 6°‘ contributions to the equations of fluid motion in the limit 


(^), when all the gradients of the metric can be neglected. Furthermore, we restrict our consideration to 
the leading in derivatives term (ISSf). Primarily, we shall be interested in the influence of the nonlinear terms 
in the equations of motion on the dynamics of the transverse modes. 

In order to find the quadratic contribution, we need to take a variation with respect to of the linearized 
equations of motion (| 6 ^ . It is convenient to express this variation in terms of 0“ rather than 5x^ just as 
it was done in deriving the linearized equations. The variations of all the structures entering the expression 


(^) are evaluated according to the rules (^), but we have to take into account that 

56^ = -e^VbO^. 


(89) 


Since the hydrodynamic equations in the Fock-Taub representation possess the property (1^), it is useful to 
develop in a Taylor series ^J\h\Ec rather than Ec- Then 


= 5{y/\h\E^c) + y\V\h\E^c) + ■■■ = + 5E^,Va9^ + ^^E^, +...), 


(90) 
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where in the last equation we suppose, for simplicity, that the background fields satisfy the classical equations 


of motion (55). The ellipses mark the higher terms of the series in 6°‘. Performing rather bulky, but simple. 


calculations, in the limit (63), we arrive at 


d 


VW\E^c =-^\ [2/3Vo(l + daen + 4/3Vo^o] (4^0 + ^c) + 2/3Vo0“<9,0a+ 


/3Vo + {daO^f + dbO^daO’^) + 2/3^^" (20o(l + daO^) + + 4/3Vo"^o1 + • • • = 0. (91) 


+ 5 


Remark that p® described in ( |63| ) with (5 = const does provide the solution to the equations of motion 
not only in the limit when the gradients of the metric can be neglected, but also for an arbitrary stationary 
metric with the Killing vector Of course, the last equality in ( | 6 ^ ) does not hold for an arbitrary stationary 
metric. The held 0“ has the mass dimension —1. It follows from (|9l|) that the amplitude of this held is small 
provided that 

|k|r ~ eyxph <c 1, (92) 

where \ph is a characteristic wavelength of phonons. To put it in another way, recollecting the dehnition 
of the deviation of huid particles from their average positions should be small in comparison with the 
wavelength. Only in this case can the higher terms of the expansion (91) be thrown away. 

Let us investigate the inhuence of the nonlinear terms on the transverse perturbations of the huid. We 
shall solve equation (91) using the perturbation theory assuming that the amplitude of the perturbation 0“ is 
small. If the transverse perturbations are small then, as we saw, their dynamics are trivial on the background 

only due to the 


(63). The small transverse modes can evolve in a nontrivial way on the background 
presence of nonlinear contributions of the longitudinal modes with a sufficiently large amplitude. At that 
the amplitude of the longitudinal modes has to satisfy 


for we can use (91) for their description. In that 
case, the longitudinal wave of a sufficiently large amplitude induces small transverse perturbations. This 
mechanism for generation of the transverse perturbations is analogous to the acoustic Rayleigh streaming 
(see, e.g., [^, 115-118|), but it does not need the viscosity of fluid to work. In fact, the vorticity (16) 


is conserved in our case. For relativistic fluids, we have to distinguish the transverse and entropy-vortex 
perturbations in higher orders of perturbation theory (see equation (98) below). 

Let A be an amplitude of the longitudinal modes 9\ and a be an amplitude of the induced perturbations 
tpa (both longitudinal and transverse). The total perturbation takes the form 


0a = el + 'ljja = 0l + 'lpl + 


(93) 


Then we deduce from (91) that 

|k|a ~ <C |k|A. (94) 

If, in the region of space at issue, the transverse perturbations are absent at the initial instant of time (that 
we shall assume) then their amplitude will be of the order (p4[) after the perturbation comes to the region 


considered. Therefore, the choice of 0a obeying (^) as the zeroth order approximation to the solution of 


(W) is justihed. 


In the leading order, equation (^) is written as 




+ 




=: 5E, 






- /c = 0, 


(95) 


Qa^'>pa 

that is, the linear in 6a terms in ( |9l|) (see also (|65D) are replaced by il^a, and the quadratic terms are 
substituted for 9\. Employing (|^, we derive from (|9l|) , 

/o = l3‘^^p'Qdij9dij9 - 2^‘^p'^di9di9 + [Mf [p'g + cl{4.j3‘^p'^ + 6p'o) + c^(4/3^p()" + 14^Vo' + Vo)] }, 
fi = ‘^P‘^Po{clA9di9 - dj9dij9), 


(96) 


where 9a = da9. Splitting the linear part of (^) as in (^), we have for the transverse part of the perturbation 
at c = i, 

= dr [dk[i9dj^k^ - ciAdii9dj]9]. (97) 
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Taking into account that the fluid flow was unperturbed at the initial instant of time r, we come to 


diiipf] = dkiiOdj^kO - c^AdiiOdj^l (98) 

This equation describes the induced transverse perturbation moving along with the sound wave. Equation 


(98) can be derived in a more straightforward manner from the vorticity potential conservation law (13). In 
the case at issue, = ta and 


ta — daOo + 9a — 29Q(da9Q + 6a) + OcdaO^ + 5j((l — 20q — (P + 40g) + . . . 




where the dots denote the terms of higher orders in Substituting (p 
derivative of ti, and setting it to zero, we come to (|98|). Actually, equation 
accurate within the second order in 

The potential part of the perturbation ij^a is specified by the system of equations 


(99) 


into (pii^), taking the external 
is the zero vorticity condition 


2/5Vo [cs ^^0 + A^/>] = /o, -2/3‘^p'qA{'iPo + ^A) = difi, 


( 100 ) 


where 'i/’i = diifj. If the fluid flow was unperturbed at the initial instant of time in the region of space 
considered then one derivative with respect to r in the second equation can be “canceled”. As a result, 
summing the equation obtained with the first one, we deduce 


2/3^Po [cs - A^/^o] = fo +difi, 


( 101 ) 


i.e., the wave equation with a source. This equation describes, in particular, the appearance of small overtones 
in a sound wave that have a doubled frequency as compared with the initial one. This phenomenon is well- 


known in acoustics (see, e.g., ||^, 117|). 

Equation ( lOlD can be considerably simplified, if one keeps there only the leading term at x —0 (see the 
asymptotics (39) and (pl|)). In this limit, 

( 102 ) 




Po 


and 


“ 5 [(if - 0 


fl 


2/3^Po 

djfi 

2/32p') 


—djOdijO, 

-di{dij6dj6). 


Whence 


(/o + ai/0/(2/3Vo) ~ A0A0 


W 


-l)=-lAeA9 = -^dr[m^] 


Consequently, 


cJ^^o-Ai;o = -i^dr[{Ae)% 


(103) 

(104) 

(105) 


These calculations show that the classical perturbation theory, i.e., in fact, the perturbation theory for 
averages, is applicable in the flat space-time limit x —)• 0 as long as the condition (92) holds. In Sect. 0, we 
shall see that the standard in-out quantum perturbation theory does not work in this limit. 

Thus, the accounting of the nonlinear terms of second order in the fluid equations of motion does not 
lead to self-maintained nontrivial dynamics of the transverse perturbations. The nontrivial dynamics of the 
transverse modes may arise at the third order. Supposing that the connected Green function (T{0“0^}) is 


not zero and neglecting the connected correlator of the third order in equations (^) expanded up to the third 
order in 9a and averaged over the vacuum state, one can derive a closed system of equations for the averaged 
field 0“ and the Green function {T{9°‘9^}). This approach is equivalent to the leading approximation in the 
effective action formalism for composite operators ||119| . In our case, the role of an independent composite 
operator is played by the two-point Green function. The study of this system of equations will be given 
elsewhere. 
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6 Concrete form of the classical potential 

In the previous sections, we have worked with a general expression ( |3^ for the potential without any 
reference to its peculiar structure and the values of the coefficient entering it. fn this section, we consider 
several particular cases of (|3^ and discuss the restrictions on the coefficients coming from the normalization 
conditions (^hj), (^Tj), and (]^. 

First of all observe that if we truncate the expansion (|3^ ) over I at the terms with Imax = L (recall we 
suppose that the anomalous scaling is small and so the coefficients at the logarithms are small as well) then 
the pressure will contain 3(L + 1) structures Hence, adjusting the coefficient in such a potential, 

we can make vanish no more than (3L + 2) derivatives with respect to at = 1. This is equivalent to a 
vanishing of (3L + 2) derivatives of the potential with respect to t‘^ at = 1. As a result, in the weak field 
limit, the polytropic index is equal to [see 




= 3L + 1. 


(106) 


fn all, the pressure (32) at Imax = L contains 6 (L + 1) arbitrary constants. Therefore, usually, it is not a 
problem even for the maximal polytropic index Umax to impose the additional normalization conditions (^^ 
and (]^). Nevertheless, having imposed the constraints ( ^7| ) and (^^, the natural electroweak scale arises 
in the potential - the normalization conditions involve the vacuum expectation value of the Higgs field rj 

So, in order to satisfy the asymptotics (t 


and the Higgs boson mass rriy-. So, in order to satisfy the asymptotics (^) with the physical value (43) of 
the constant A at small n < 8 , we have to fine-tune the constants bkis in such a way that the coefficients 
at In^ agree with the value (^^ at (p = rj. We shall see this below for the concrete examples with 
L = {0,1,4}. 

Let us start with the simplest model L = 0. fn this case, Umax = 1 and coincides with the minimal value 
of n admissible by the conditions 1 and 2 of Sect. |^. Imposing the normalization conditions (36), (^), and 


we arrive at 


Po = 


m; 

8 ry 


r? 


2\2 


5100 


- 1)^ 


(107) 


i.e., the potential falls into two noninteracting parts - the Higgs potential and the potential (the pressure) 
of the field The value of the constant 6 ioo is fixed by the condition 


bioo = — Al¬ 


im 


The value of Ai is given in (^4|). 

Now we consider L = 1. The model of this type with the additional quantum corrections was studied 
in 1^. fn this case, Umax = 4 and we put n = Umax- Imposing the normalization conditions providing the 
polytropic index n = 4 in ( ^^ and the normalization conditions (^), ( |38|) (eight conditions, in all), we 
obtain a four-parameter family of potentials 


Po = 


ll^|2 .„2^2 I t\j,\2 .„2^ [^010^^2 i .l2^ , ^HO 


3ri 




— &210 


- l)i3t^ - 2 + - {t^ + + (109) 


In particular, for \(p\ = rj, we have 


^210 , 


= &21o[ 3 - 3f^ + + 4f2 + 1) lnf2] = _^(f2 _ ^ 


( 110 ) 


The constant 6210 is fixed by the asymptotics (3£) with the constant A presented in (44). 

As the last example, we consider the case L = 4. fn 0 , it was suggested that we can get rid of 
the dimensional constant in the Higgs potential (see the notation in |^) replacing this constant by at"^. 
Then, redefining t —>■ only two dimensional constants remain in the standard model with gravity: the 

gravitational constant and the cosmological one. Let us show that such a model (for definiteness, we call it 
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conformal) can be realized at L > 4. Namely, in the conformal model, the pressure 
the dimensional constants, after the dilatation t 


and 




4-2fc 


does not contain 


( 111 ) 


In all, there are 3{L + 1) arbitrary constants in the potential, the coefficient at being of the order 

we need to adjust the 


of the electroweak scale. Therefore, in order to arrive at the asymptotics (3£), 

coefficients in the potential such that the polytropic index n > 8 . Taking into account the requirements (| 
(38), we have (n + 4) > 12 conditions. For L = 3, the pressure contains 12 constants, but the system proves 
to be inconsistent. Actually, this system of conditions does not possess any solution already at n > 7. Only 
at L > 4 can we realize the conformal model. 

Let us consider L = 4. It turns out that the natural polytropic index is n = 8 for this model. For n >9, 
having imposed the normalization conditions, the coefficients at the logarithmic corrections in (^) prove to 
be large in comparison with the same structures without logarithms. This contradicts to the requirement 
of a small anomalous scaling that we assume. For n = 8 , the normalization conditions (^^, (|3^ , imply a 
three-parameter family of potentials 


Po = 


9 2 

m^?7 ^ 


(r^ - - lOtV^ -P 3r^) Int^ + M79t^ -P IGtV^ - 3l/) In^ t^-P 


-P 4(13t^ - 88r^)t^ In^ “ ie + ^03^"^ 24(r^ - Inf + 9{f - f) In^ f + (r'^ - Sff + f) In^ f 




+ 48(r^ -fflnf + l + IQff - 5lf) In^ f - | ( 8 r^ - f)f In^ f - - ff) In^ f 


+ ^24 


168(r" - ff In f + - t") (48 In" f - In^ f) - 72ff In^ f 


+ 

( 112 ) 


where r := \4>\/rj. The asymptotics 
constants. For instance, 


^03 = 


with the constant A given in (44) fixes one of the three arbitrary 
35 7m" 36 i4 15524 


57122 384r/2 




rj-i 


(113) 


and then 


Po = ■ 

_1_ ^ 5m| 


2 2 
m^rj 


315 


5127?'^ 57122 

2 


^57122 384r;: 

'in 4 n.2 2 


- 5147 /" 24(r" - t")" Inf" + - 6 t"r" - ff) In" f + (3r^ - 6 t"r" + ft") In" + (^ - ff) In^ t 


28561. 
2+2 I /Q 4 


57122 512??- 

2 

X 

768n 


53m| 4 

- J 


f *")] + [( 

315 ^,^4^ , 2 .2 , r,' 35 _ _ 140 n,2 2 , / 35 _ w41 , 3 n _ 

A jm 7 “T LI 57122 3Mn^A “T Vs^ 2856 iF ' ■rV57l22 768r)^ A J '' 


)r- - 


3+4i 


768??' 

+ 2 „ 2 i 


■In^ f[- 


A+2 


— 524 


192(r" - fy Inf + {f - t^)(871n" f + In^ f) + (15r^ - 48t"r" + 15t^) In" f . (114) 


One can check that the coefficients at the logarithms are small for reasonable, in particular, for vanishing, 
values of the constants 5 i 2 , 524 . Since the polytropic index is an even number, the domains with a negative 
energy density do not appear in the weak field limit. 

The question as to make the gravitational sector of a theory invariant with respect to the global conformal 
(Weyl) transformations by substituting the Hilbert-Einstein action for the terms with second derivatives ( |^ 
and the logarithmic corrections to them remains open. As we have already noted in Sect. we may hope this 
theory is conformally invariant on the quantum level. This happens when the coefficients at the logarithms 
in (nd) coincide with the renormalization group functions (3 and 7 taken on the cutoff scale A. Of course, as 
the expression (114) is approximate, the coincidence should be approximate too. The analysis above shows 
that the conformal model of a dark matter can be realized for the polytropic indices n > 8 only. Therefore, 
it can be excluded by the experimental data regarding the properties of a dark matter. 
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7 Some quantum effects 


In this section, we try to estimate the decay rate of the Higgs boson to phonons and evaluate the one-loop 
contribution of phonons to the effective potential. 

To begin with, we consider the one-loop correction to the effective action coming from the quantum 
fluctuations of phonons on the stationary background with the Killing vector in the leading order in 
derivatives. If the fluid rests in the adapted system of coordinates, where = (1,0,0,0), then the phonon 
dispersion law has the form 

+ 2csg^^u)ki -|- clg'^^kikj = 0. (115) 

In the flat space-time limit, 

ijQ,<A 


= -- 

2 




dk 


2 ( 27 r) 


r0(A — Cs|k|)cs|k| = — J d^xc 


„-3 


IOtt^ 


(116) 


When the derivatives of the metric can be neglected, the dependence on it is easily recovered (for details, 
see (18) in p) 


p(i) = _ 


J d^x^/\^\c 


-3 


IGvr^ 


(117) 


The explicit expression for Cg is given by formula (|4lD . The expression obtained is singular in the flat space- 
time limit when Cg —?• 0. In order to comply with the conditions 1 and 2 of Sect. this contribution to 
the effective action must be completely canceled out. One can think that it is possible to cancel only the 


singular in the flat space-time limit part of the contribution (07). This can be done, if one adds to ( 117 ) 
the terms proportional to 

(t^ - l)-3/^ (t^ - 1)-V2^ (t2 _ 1)1/2^ 

with the corresponding coefficients. However, upon stretching t —>■ A~^t and developing as a series in 
these terms lead to the appearance of the coupling constants with arbitrary large negative mass dimension 
in the theory. This contradicts the property mentioned in Sect. after the condition 3 that we wish to 
preserve. 

In conclusion, we consider the decay of the Higgs boson to phonons. As long as the pressure ( p^ contains 
the terms that depend on the product of the fields (j) and one may suppose that this decay is possible. 
Then the experimental data concerning the Higgs boson decay can be used to impose the constraints on the 
values of the coefficients at the corresponding terms. However, we shall show now that, in the weak field 
limit, the standard in-out perturbation theory is not applicable to phonons since they constitute a condensate 
- the macroscopic field - and, in a certain sense, are not elementary fields of the theory in this limit. 

Consider the case when the gradients of the metric can be neglected and the background field obeys ( |^ . 
We shall consider the longitudinal phonon modes only. It follows from (|65|) that the canonically normalized 
field is defined as 

»o“:=(-^)‘'> = (|)‘''>. (119) 

The expressions for variations of the potential (^^ with respect to 0“ are present in (^) and (^). It is 
only relevant for us that, upon variation, the field 0“ always enters with one derivative. Let us show that 
the Higgs boson decay to a larger number of phonons is more probable than to a lesser number of them. In 
other words, the decay process is not perturbative. 

Indeed, according to the standard expression (see, e. 
particles in the rest frame of a decaying particle, we have 

\M\'^_/ci\n^ dk 


g., |120|) for the differential decay rate to n identical 


dT{x nO) = (27r)^(5^ fp — ^ k. 


2=1 


2mvra! 


n 


2a;i(27r)3’ 
2=1 ^ ' 


( 120 ) 


where = (m^,0,0,0) is the four-momentum of the Higgs boson, oji = Cs|kj|, and M is the invariant 
scattering amplitude. Its square is of the order 


|A/|^ 


V 


’/ d 

W 9po ’ 

V5p2^ 

03 

to 




( 121 ) 


'l=»?.p =1 i=i 
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where Wj come from the derivatives with respect to r acting on the held 9^. To obtain the total decay rate 
for X n0, it is necessary to integrate (|120| ) over all the phonon momenta. Roughly, this integral can be 
estimated as follows: every k* and dkj are replaced by and the delta-functions produce the factor 


Then 


r(x ^ nO) 


m|n! 


d 

dp^ 


^ dpo 


>|=r),p2=l 


m" 






( 122 ) 


There also exits the contribution of the same order from the variation of Y^[h[ and expansion of the Higgs 
held in 9°‘. Of course, the dependence on n in this formula is incorrect since we did not evaluate the 

integrals exactly. Nevertheless, this formula implies that for 


ml 




> 1 , 


(123) 


the standard perturbation theory is not applicable. Taking into account the actual value of £q in the weak 
held limit (d2), we see that to handle the decay process x by means of the standard S-matrix technique 


is a nonsense. The condition (1^) agrees with the general strong coupling condition found in 

In many respects, this situation resembles the theory of a quantum held with the held strength renormal¬ 
ization constant Z —>• 0 (see, for details, |121|). The kinetic term disappears from the action at Z —?■ 0 for 
that held, just as in our model of phonons. In this sense the phonons are not elementary helds of the theory 
in the weak held limit. Only in the strong gravitational helds do the quantum properties of these particles 
become relevant. 

On the other hand, we saw in Sect. |^that the classical perturbation theory, i.e., the leading order of the 
perturbation theory for averages, works even in the weak held limit provided (|^ ) is fulhlled. This allows us 
to hope that the properly formulated in-in perturbation theory can be applied to describe quantum dynamics 
of the held In particular, it follows from (|6^ and (^) that, in the leading order, the wave packet of 
Higgs bosons inhuences the held 0“ solely though the sound speed Cs due to small variations of the vacuum 
expectation value p. So, we have for the perturbation of the phonon potential 56 over the background 0, 


- A59 = 


- 2 \ 


(124) 


Consequently, if the phonons are absent at the initial instant of time, they will not be produced later on. In 
fact, this means that the Higgs boson and the phonons are not mixed. 


8 Conclusion 

So, we have constructed a self-consistent formalism of quantum held theory on a curved background that 
solves dynamically the problem of dependence of observables on a choice of the vacuum state for quantum 
helds. At the same time, we have to state that many hopes formulated in regarding the predictive power 
of the theory developed are not justihed because of the perturbative nonrenormalizability of relativistic 
hydrodynamics and its non-perturbativity at large momenta. A virtual lack of reliable experimental data 
that could reveal the additional peculiarities (symmetries) of the huid effective action leaves a large freedom 
in a selection of the concrete model. However, we saw that, in the weak gravitational held limit, one can 
establish certain properties of this model that allow one to falsify it so long as this huid is identihed with 
a considerable part of a dark matter. Moreover, we can use the derivative expansion and construct the 


standard perturbation theory for processes with the energies satisfying the inequality inverse to ( 123 ), where 
ought to be replaced by the energy of a process. In particular, it would be interesting to consider the 
decay of a low energy graviton to phonons when this process is perturbative. This process describes the 
inhuence of the dark matter on the graviton dynamics and seems result in a hnite lifetime of the graviton. 
The classical perturbation theory is applicable under the much weaker condition (p^). 

In order to identify the quantum gravitational anomaly with a dark matter or its considerable part, 
further investigations regarding the refinement and interpretation of the astrophysical data are also needed. 
This will allows us to find the universal constants A and n characterizing the fluid equation of motion in 
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the limit of a weak gravitational field, or to prove that the dark matter is not described in this limit by the 
universal polytropic equation of state. The latter will signify that the hypothesis about the identification of 
the quantum gravitational anomaly with the dark matter is not valid. This, however, will not mean that 
the gravitational anomaly is absent. The relativistic fluid described by the field arises with a necessity 
in quantum field theory with gravity provided that the basic principles of quantum field theory and general 
relativity hold and the general covariance (the background independence) remains intact on a quantum level, 
i.e., the general covariance is a fundamental symmetry of Nature. As we have already mentioned, this is 
the most natural solution to the problem of time in quantum gravity or, put another way, the problem of a 
unique representation of the algebra of observables (other similar models see in p0|-22, p^). If something 
can be called the flow of a physical (not coordinate) time then the vector field would be the most natural 
candidate to this role. 
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A Linearized equations of motion 

In this appendix, we linearize the expressions ( p^ with respect to the small perturbation In fact, as we 
have already mentioned in the main text, we need to evaluate the Lie derivative of the expressions ( |54[ ) using 
formulas (^3|). The outcome of these calculations is as follows 

5El = 2Va[t^t^R{^(A) - 4tc/Vb0d) + + 2tcERVabe’’ - O^VaR^ct^ + 2RVc{tH^)Vaeb, 

5El = AVd{t^th%e^VsRab + 2R,i,Vae^)pc + RabP^i^icOs) - QPcg^^dOs)] ] - 2Vb[tH\e^VsRac+ 

+ RscVaO^ + RasVcO^ " ARabt<^P^VdOs)] " sRab + 2Rsb^a9^ " Vrf^s)V,(t“t'>)- 

- 2Racet^Vbse^ + 4R„fet“t^(Ve(t'^/)Vrf0, + t^Pc^dsO^), 

[vH^V^aOb) - 2V\tbEp^Va9s) - V“(V(,0,)V*4) + 

+ V^E{9^Rdisba) + ^isa)9b) + V^(V,rf0fe + 9^Rasbd)g^] - 2V\Vc{th^p^)Vs9d- 

- 2evas9H^VHbPc + eVas9^V\ + Va (V(,0fe)- 2t^p^Vb9sV\ - 2t^V\V^,9s)P^- 

- V*(V(,0b)V''4) - 2V\Efp^Vs9b) + Vbs9^vHc + VH^{9^Rb(dcs) + ^(ds)9c)+ 

+ V\Vbs9c + 9<^Rdbcsy) - ^et^Pc [^sdG^^VHb - V^(V(,0rf)V‘='4) - 2V\ht^p^Vs9d)- 

- At<^p^Vd9sS/\ + ^ib9s)^^t^ + VH^{9^Rridbs) + ^ids)9b) + ^^{^sd9b + 9^Rrsbd)t^ 

^ 4{v,4it“fVfe,r - EVb{t’^Ep'^V,9d)] - V,{efp^Vabt^)^s9d - 


6Ej = 2\ Vet 


6E^ = 


- ApeVabth^^P^^dSr " 2peVabit^^^p'^^d9r) + Pc^bd9^))] - Pc^abt^^ds9^] , 
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5El = -2{Vet[V2(tVVa0fe) + V^(V(,0fe)V“t) - V^tVabO^] + V2t[Ve(tV)V„0fe + t“neVab0^] + 

+ Va[vh{Ev!>v^A) - 2,Euct<^p^Vdes) - t^uc{v\t<^p^Vdes) + - v'^tv^.r)]}, 

{Ve[a;at“tVfe(t‘='/5^Vrf0,) - 2(u;at“)2tVVd0.] - {t^Pc^dsO^ - t<^p^V,des)[yb{t^0Jan + 2{u:at^f]- 
t'^P^i^icOs) - 6fpcVs9r){Vb{t^0Jan + 2(w,t“)2) + p.t’^Vb,9^0jX + Vfe [tVVa(t''p^Vrf0,)] ] } , 
{Vp'' [Vab{t^tH^Vd9s) + VHHVab9d + 9^Rsbda)] } + 4V,fe [p<^p’^VciM^Vd9s)] - 

- Va.rVb(pvVet") + Vb{p^p^Vdt^)iVca9^ + - V4/5VVet"Vfe,r+ 

+ //5Vei^(VMr - R\ds9n - p‘^p^Vdt\Vci,9‘^ + R\sc 91] + 4.p‘^p^Vr{thH‘^Vd9s)R\cb- 

- p^p^Vdt^CeR\,i, + AVd{thH^V^,9,)Vab{p^P^) - 6t\ft^Vd9sVab{p'^p’’) + Mc[Vas9^Vb{p‘^p^)+ 


5El = -2 


- V, 


5El = V, 


+ - R\sr91)\ \ + AtH%Vds9^yab{p"p^) " UWVcd9ab{p^P^) 


2+d4 


4.2+dj.s\ 


5El = -4{Ve[pVfe(/Va(t2t‘'t"VA))] + V4tVt^Va(/VfeP^)(V(,0,) -6tVeV,0,) + 

+ tH%Va{p^P^Vl,s9^) + t\t^Vas9^VbP^] + tH^{EVds9^ - t"Verf0.)Va(p“Vb4)}. 

(125) 

All the calculations were carried out by hand and so the mistakes are possible, but unlikely. As expected, 
the first two expressions vanish at i? = 0 and Rab = 0, respectively. 
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